Discovery and characterization of novel DNA viruses in Apis mellifera: expanding the honey bee virome through metagenomic analysis.

To date, many viruses have been discovered to infect honey bees. In this study, we used high-throughput sequencing to expand the known virome of the honey bee, Apis mellifera, by identifying several novel DNA viruses. While the majority of previously identified bee viruses are RNA, our study reveals nine new genomes from the Parvoviridae family, tentatively named Bee densoviruses 1 to 9. In addition, we characterized a large DNA virus, Apis mellifera filamentous-like virus (AmFLV), which shares limited protein identities with the known Apis mellifera filamentous virus. The complete sequence of AmFLV, obtained by a combination of laboratory techniques and bioinformatics, spans 152,678 bp. Linear dsDNA genome encodes for 112 proteins, of which 49 are annotated. Another large virus we discovered is Apis mellifera nudivirus, which belongs to a group of Alphanudivirus. The virus has a length of 129,467 bp and a circular dsDNA genome, and has 106 protein encoding genes. The virus contains most of the core genes of the family Nudiviridae. This research demonstrates the effectiveness of viral binning in identifying viruses in honey bee virology, showcasing its initial application in this field.IMPORTANCEHoney bees contribute significantly to food security by providing pollination services. Understanding the virome of honey bees is crucial for the health and conservation of bee populations and also for the stability of the ecosystems and economies for which they are indispensable. This study unveils previously unknown DNA viruses in the honey bee virome, expanding our knowledge of potential threats to bee health. The use of the viral binning approach we employed in this study offers a promising method to uncovering and understanding the vast viral diversity in these essential pollinators.

A new member of the Nudiviridae from the Florida stone crab (Menippe mercenaria).

Menippe mercenaria, the Florida stone crab, supports an unconventional fishery across the southern USA and Caribbean that involves claw-removal and the return of de-clawed animals to the sea. We provide pathological, ultrastructural, and genomic detail for a novel hepatopancreatic, nucleus-specific virus - Menippe mercenaria nudivirus (MmNV) - isolated from M. mercenaria, captured during fisheries-independent monitoring. The virus has a genome of 99,336 bp and encodes 84 predicted protein coding genes and shows greatest similarity to Aratus pisonii nudivirus (ApNV) (<60% protein similarity and 31 shared genes of greatest similarity), collected from the Florida Keys, USA. MmNV is a member of the Gammanudivirus genus (Naldaviricetes: Lefavirales: Nudiviridae). Comparisons of virus genome size, preferred host environment, and gene number revealed no clear associations between the viral traits and phylogenetic position. Evolution of the virus alongside the diversification of host taxa, with the potential for host-switching, remain more likely evolutionary pathways.

Variants in the mitochondrial genome sequence of Oryctes rhinoceros (Coleoptera: Scarabaeidae) infected with Oryctes rhinoceros nudivirus in oil palm and coconut plantations.

The CRB (coconut rhinoceros beetle) haplotype was classified into CRB-S and CRB-G, based on the presence of single nucleotide polymorphisms (SNPs) in the mitochondrial cox1 gene. Mitochondrial genomes (mitogenomes) are the most widely used genetic resources for molecular evolution, phylogenetics, and population genetics in relation to insects. This study presents the mitogenome CRB-G and CRB-S which were collected in Johor, Malaysia. The mitogenome of CRB-G collected from oil palm plantations in 2020 and 2021, and wild coconut palms in 2021 was 15,315 bp, 15,475 bp, and 17,275 bp, respectively. The CRB-S was discovered in coconut and oil palms in 2021, and its mitogenome was 15,484 bp and 17,142 bp, respectively. All the mitogenomes have 37 genes with more than 99% nucleotide sequence homology, except the CRB-G haplotype collected from oil palm in 2021 with 89.24% nucleotide sequence homology. The mitogenome of Johor CRBs was variable in the natural population due to its elevated mutation rate. Substitutions and indels in cox1, cox2, nad2 and atp6 genes were able to distinguish the Johor CRBs into two haplotypes. The mitogenome data generated in the present study may provide baseline information to study the infection and relationship between the two haplotypes of Johor CRB and OrNV in the field. This study is the first report on the mitogenomes of mixed haplotypes of CRB in the field.

Atomic structure of a nudivirus occlusion body protein determined from a 70-year-old crystal sample.

Infectious protein crystals are an essential part of the viral lifecycle for double-stranded DNA Baculoviridae and double-stranded RNA cypoviruses. These viral protein crystals, termed occlusion bodies or polyhedra, are dense protein assemblies that form a crystalline array, encasing newly formed virions. Here, using X-ray crystallography we determine the structure of a polyhedrin from Nudiviridae. This double-stranded DNA virus family is a sister-group to the baculoviruses, whose members were thought to lack occlusion bodies. The 70-year-old sample contains a well-ordered lattice formed by a predominantly α-helical building block that assembles into a dense, highly interconnected protein crystal. The lattice is maintained by extensive hydrophobic and electrostatic interactions, disulfide bonds, and domain switching. The resulting lattice is resistant to most environmental stresses. Comparison of this structure to baculovirus or cypovirus polyhedra shows a distinct protein structure, crystal space group, and unit cell dimensions, however, all polyhedra utilise common principles of occlusion body assembly.

Complete genome sequence of an oryctes rhinoceros nudivirus isolated from Korean rhinoceros beetles (Trypoxylus dichotomus) in Korea.

Nudivirus-infected Korean rhinoceros beetles (Trypoxylus dichotomus) were first identified in 2015, and while a complete genome sequence of the virus has long been uploaded to the NCBI database, it has not been examined in detail. Here, we describe the genomic characteristics of Trypoxylus dichotomus nudivirus (TdNV), which represents a new Oryctes rhinoceros nudivirus (OrNV) strain, isolated from infected T. dichotomus in the Republic of Korea. We examined factors derived by the cross-species infection of OrNV from nucleotide levels to the whole genome level. Our genomic analysis study suggests that TdNV-Korea is highly conserved with other OrNVs in terms of genomic structures and genome size. Our investigation of the genomic structure revealed that TdNV-Korea has the least number of open reading frames (ORFs) of all available OrNV genomes; three hypothetical genes were notably absent only in TdNV-Korea. In addition, the genomic alteration of the nudivirus core genes discloses that various amino acid mutations caused by single-nucleotide polymorphism and short indels (insertion/deletion) were found in most of the nudivirus core genes of TdNV-Korea. Our findings provide a valuable resource for those seeking a greater understanding of cross-species nudivirus transmission and will certainly provide valuable insight for reconstruction and reinterpretation of future and previously identified OrNV strains.

Rediscovering "Baculovirus-A" (Johnson, 1976): The complete genome of 'Callinectes sapidus nudivirus'.

Callinectes sapidus, or the 'blue crab', supports an extensive east-coast USA fishery and was one of the first crustacean species in which viruses were observed. Pioneering research by Dr Phyllis Johnson led to these initial discoveries, one of which included the discovery of a virus termed "Baculovirus-A". This virus was considered a potential member of the Baculoviridae, Nimaviridae, or Nudiviridae, in which all viral members are rod-shaped dsDNA viruses found in the nucleus of their host cell. With the availability of genomic and bioinformatic tools, such as Illumina HiSeq and assembly programs, it is now possible to assemble the genomes of viruses and gain additional genomic insight, which can shed light on viral taxonomy. Using these tools, alongside electron micrographs and histology slides, we reveal that the hepatopancreas-infecting 'Baculovirus-A' from Callinectes sapidus is a member of the Nudiviridae, resembling genetic and protein similarity to other crab and lobster infecting nudiviruses from the Gammanudivirus genus. Histologically, the virus causes nuclear hypertrophy as observed for other gammanuriviruses. The genome of the virus is circular, 122,436 bp in length, and encodes a predicted 98 protein coding genes, including all of the nudivirus core genes. The prevalence of virus from across Florida, USA, is provided alongside a genomic comparison of the new viral genome against other Gammanudivirus species, revealing the average prevalence to be 2.2% and that Callinectes sapidus nudivirus is distantly similar to the recently described Carcinus maenas nudivirus from Canada.

Conserved Viral Transcription Plays a Key Role in Virus-Like Particle Production of the Parasitoid Wasp Venturia canescens.

Nudiviruses are large double-stranded DNA viruses related to baculoviruses known to be endogenized in the genomes of certain parasitic wasp species. These wasp-virus associations allow the production of viral particles or virus-like particles that ensure wasp parasitism success within lepidopteran hosts. Venturia canescens is an ichneumonid wasp belonging to the Campopleginae subfamily that has endogenized nudivirus genes belonging to the Alphanudivirus genus to produce "virus-like particles" (Venturia canescens virus-like particles [VcVLPs]), which package proteic virulence factors. The main aim of this study was to determine whether alphanudivirus gene functions have been conserved following endogenization. The expression dynamics of alphanudivirus genes was monitored by a high throughput transcriptional approach, and the functional role of lef-4 and lef-8 genes predicted to encode viral RNA polymerase components was investigated by RNA interference. As described for baculovirus infections and for endogenized nudivirus genes in braconid wasp species producing bracoviruses, a transcriptional cascade involving early and late expressed alphanudivirus genes could be observed. The expression of lef-4 and lef-8 was also shown to be required for the expression of alphanudivirus late genes allowing correct particle formation. Together with previous literature, the results show that endogenization of nudiviruses in parasitoid wasps has repeatedly led to the conservation of the viral RNA polymerase function, allowing the production of viruses or viral-like particles that differ in composition but enable wasp parasitic success. IMPORTANCE This study shows that endogenization of a nudivirus genome in a Campopleginae parasitoid wasp has led to the conservation, as for endogenized nudiviruses in braconid parasitoid wasps, of the viral RNA polymerase function, required for the transcription of genes encoding viral particles involved in wasp parasitism success. We also showed for the first time that RNA interference (RNAi) can be successfully used to downregulate gene expression in this species, a model in behavioral ecology. This opens the opportunity to investigate the function of genes involved in other traits important for parasitism success, such as reproductive strategies and host choice. Fundamental data acquired on gene function in Venturia canescens are likely to be transferable to other parasitoid wasp species used in biological control programs. This study also renders possible the investigation of other nudivirus gene functions, for which little data are available.

A high-quality de novo genome assembly based on nanopore sequencing of a wild-caught coconut rhinoceros beetle (Oryctes rhinoceros).

BACKGROUND: An optimal starting point for relating genome function to organismal biology is a high-quality nuclear genome assembly, and long-read sequencing is revolutionizing the production of this genomic resource in insects. Despite this, nuclear genome assemblies have been under-represented for agricultural insect pests, particularly from the order Coleoptera. Here we present a de novo genome assembly and structural annotation for the coconut rhinoceros beetle, Oryctes rhinoceros (Coleoptera: Scarabaeidae), based on Oxford Nanopore Technologies (ONT) long-read data generated from a wild-caught female, as well as the assembly process that also led to the recovery of the complete circular genome assemblies of the beetle's mitochondrial genome and that of the biocontrol agent, Oryctes rhinoceros nudivirus (OrNV). As an invasive pest of palm trees, O. rhinoceros is undergoing an expansion in its range across the Pacific Islands, requiring new approaches to management that may include strategies facilitated by genome assembly and annotation. RESULTS: High-quality DNA isolated from an adult female was used to create four ONT libraries that were sequenced using four MinION flow cells, producing a total of 27.2 Gb of high-quality long-read sequences. We employed an iterative assembly process and polishing with one lane of high-accuracy Illumina reads, obtaining a final size of the assembly of 377.36 Mb that had high contiguity (fragment N50 length = 12 Mb) and accuracy, as evidenced by the exceptionally high completeness of the benchmarked set of conserved single-copy orthologous genes (BUSCO completeness = 99.1%). These quality metrics place our assembly ahead of the published Coleopteran genomes, including that of an insect model, the red flour beetle (Tribolium castaneum). The structural annotation of the nuclear genome assembly contained a highly-accurate set of 16,371 protein-coding genes, with only 2.8% missing BUSCOs, and the expected number of non-coding RNAs. The number and structure of paralogous genes in a gene family like Sigma GST is lower than in another scarab beetle (Onthophagus taurus), but higher than in the red flour beetle (Tribolium castaneum), which suggests expansion of this GST class in Scarabaeidae. The quality of our gene models was also confirmed with the correct placement of O. rhinoceros among other members of the rhinoceros beetles (subfamily Dynastinae) in a phylogeny based on the sequences of 95 protein-coding genes in 373 beetle species from all major lineages of Coleoptera. Finally, we provide a list of 30 candidate dsRNA targets whose orthologs have been experimentally validated as highly effective targets for RNAi-based control of several beetles. CONCLUSIONS: The genomic resources produced in this study form a foundation for further functional genetic research and management programs that may inform the control and surveillance of O. rhinoceros populations, and we demonstrate the efficacy of de novo genome assembly using long-read ONT data from a single field-caught insect.

The Complete Genome of Chelonus insularis Reveals Dynamic Arrangement of Genome Components in Parasitoid Wasps That Produce Bracoviruses.

Bracoviruses (BVs) are endogenized nudiviruses in parasitoid wasps of the microgastroid complex (family Braconidae). Microgastroid wasps have coopted nudivirus genes to produce replication-defective virions that females use to transfer virulence genes to parasitized hosts. The microgastroid complex further consists of six subfamilies and ∼50,000 species but current understanding of BV gene inventories and organization primarily derives from analysis of two wasp species in the subfamily Microgastrinae (Microplitis demolitor and Cotesia congregata) that produce M. demolitor BV (MdBV) and C. congregata BV (CcBV). Notably, several genomic features of MdBV and CcBV remain conserved since divergence of M. demolitor and C. congregata ∼53 million years ago (MYA). However, it is unknown whether these conserved traits more broadly reflect BV evolution, because no complete genomes exist for any microgastroid wasps outside the Microgastrinae. In this regard, the subfamily Cheloninae is of greatest interest because it diverged earliest from the Microgastrinae (∼85 MYA) after endogenization of the nudivirus ancestor. Here, we present the complete genome of Chelonus insularis, which is an egg-larval parasitoid in the Cheloninae that produces C. insularis BV (CinsBV). We report that the inventory of nudivirus genes in C. insularis is conserved but are dissimilarly organized compared to M. demolitor and C. congregata. Reciprocally, CinsBV proviral segments share organizational features with MdBV and CcBV but virulence gene inventories exhibit almost no overlap. Altogether, our results point to the functional importance of a conserved inventory of nudivirus genes and a dynamic set of virulence genes for the successful parasitism of hosts. Our results also suggest organizational features previously identified in MdBV and CcBV are likely not essential for BV virion formation. IMPORTANCE Bracoviruses are a remarkable example of virus endogenization, because large sets of genes from a nudivirus ancestor continue to produce virions that thousands of wasp species rely upon to parasitize hosts. Understanding how these genes interact and have been coopted by wasps for novel functions is of broad interest in the study of virus evolution. This work characterizes bracovirus genome components in the parasitoid wasp Chelonus insularis, which together with existing wasp genomes captures a large portion of the diversity among wasp species that produce bracoviruses. Results provide new information about how bracovirus genome components are organized in different wasps while also providing additional insights on key features required for function.

The naked truth: An updated review on nudiviruses and their relationship to bracoviruses and baculoviruses.

Nudiviruses (Nudiviridae) are double-stranded DNA viruses with enveloped and rod-shaped virions. Several insect orders (e.g., Diptera, Lepidoptera, Coleoptera, Orthoptera) and aquatic crustaceans are susceptible to nudivirus infections, which can result in varied degrees of disease in all developmental host stages. Their pathogenicity endangers insect rearing and crustacean aquacultures, but has also proven effective in biocontrol against Oryctes rhinoceros infestations. This literature review aims to present all known nudivirus species and provide a comprehensive Nudiviridae phylogeny by including recently described nudiviral isolates, and discuss this phylogeny in comparison to current opinions and taxonomical propositions. Moreover, we aim to clarify biological, pathological and genomic differences or similarities between nudiviruses and related entomopathogenic viruses, including baculoviruses (Baculoviridae) and bracoviruses (Polydnaviridae). A phylogenetic analysis using 17 concatenated nudivirus core genes resulted in the expected structure with the genera Alphanudivirus and Betanudivirus, as well as the most recently recognized genera Gammanudivirus and Deltanudivirus. The hymenopteran Osmia cornuta nudivirus (OcNV) groups closest with the hymenopteran Fopius arisanus endogenous nudivirus (FaENV) and does not share a most common ancestor with the hymenopteran bracoviruses. Except for one node, all clades are highly supported. The proposition of a recent study to assign subgroups to the alphanudiviruses might be legitimate, but more hymenopteran and orthopteran nudiviruses, especially in bees and cricket, need to be identified to resolve this proposal. In addition, freshwater and marine nudiviruses might form taxonomic subgroups among gammanudiviruses as well, but more aquatic nudiviruses need to be identified and sequenced for better resolution. Furthermore, the search for nudiviruses in insects with (semi)aquatic life stages may aid in finding the missing link that led to the manifestation of aquatic nudiviruses.

Diverse Host Immune Responses of Different Geographical Populations of the Coconut Rhinoceros Beetle to Oryctes Rhinoceros Nudivirus (OrNV) Infection.

Incursions of the coconut rhinoceros beetle (CRB), Oryctes rhinoceros, into different islands in the South Pacific have been detected in recent years. It has been suggested that this range expansion is related to an O. rhinoceros haplotype reported to show reduced susceptibility to the well-established classical biocontrol agent, Oryctes rhinoceros nudivirus (OrNV). Our understanding of the genetic characteristics which distinguish the population of O. rhinoceros that has recently established in Solomon Islands from other well-established populations across the region is very limited. Here, we hypothesized that the recently established O. rhinoceros population should have greater innate immune responses when challenged by OrNV than those of well-established and native O. rhinoceros populations. We used the RNA sequencing (RNA-Seq) approach to generate gene expression profiles of midgut tissue from OrNV-infected and noninfected individuals collected in the Solomon Islands (recent incursion), Papua New Guinea and Fiji (previously established), and the Philippines (within the native range). The collections included individuals from each of the three major mitochondrial lineages (CRB-G, CRB-PNG, and CRB-S) known to the region, allowing us to explore the specific responses of each haplotype to infection. Although insects from the Philippines and Solomon Islands that were tested belong to the same mitochondrial lineage (CRB-G), their overall responses to infection were different. The number of differentially expressed genes between OrNV-infected and noninfected wild-caught individuals from the four different locations varied from 148 to 252. Persistent OrNV infection caused a high level of induced antimicrobial activity and immune responses in O. rhinoceros, but the direction and magnitude of the responses were population specific. The insects tested from the Solomon Islands displayed extremely high expression of genes which are known to be involved in immune responses (e.g. coleoptericin, cecropin, and serpin). These variations in the host immune system among insects from different geographical regions might be driven by variations in the virulence of OrNV isolates, and this requires further investigation. Overall, our current findings support the importance of immunity in insect pest incursion and an expansion of the pest's geographic range. IMPORTANCE Oryctes rhinoceros nudivirus (OrNV) is a double-stranded DNA (dsDNA) virus which has been used as a biocontrol agent to suppress coconut rhinoceros beetle (CRB) in the Pacific Islands. Recently a new wave of CRB incursions in Oceania is thought to be related to the presence of low-virulence isolates of OrNV or virus-tolerant haplotypes of beetles (CRB-G). Our comparative analysis of OrNV-infected and noninfected CRBs revealed that specific sets of genes were induced by viral infection in the beetles. This induction was much stronger in beetles collected from the Solomon Islands, a newly invaded country, than in individuals collected from within the beetle's native range (the Philippines) or from longer-established populations in its exotic range (Fiji and Papua New Guinea [PNG]). Beetles from the Philippines and the Solomon Islands that were tested in this study all belonged to the CRB-G haplotype, but the country-specific responses of the beetles to OrNV infection were different.

Identification and Full Characterisation of Two Novel Crustacean Infecting Members of the Family Nudiviridae Provides Support for Two Subfamilies.

Multiple enveloped viruses with rod-shaped nucleocapsids have been described, infecting the epithelial cell nuclei within the hepatopancreas tubules of crustaceans. These bacilliform viruses share the ultrastructural characteristics of nudiviruses, a specific clade of viruses infecting arthropods. Using histology, electron microscopy and high throughput sequencing, we characterise two further bacilliform viruses from aquatic hosts, the brown shrimp (Crangon crangon) and the European shore crab (Carcinus maenas). We assembled the full double stranded, circular DNA genome sequences of these viruses (~113 and 132 kbp, respectively). Comparative genomics and phylogenetic analyses confirm that both belong within the family Nudiviridae but in separate clades representing nudiviruses found in freshwater and marine environments. We show that the three thymidine kinase (tk) genes present in all sequenced nudivirus genomes, thus far, were absent in the Crangon crangon nudivirus, suggesting there are twenty-eight core genes shared by all nudiviruses. Furthermore, the phylogenetic data no longer support the subdivision of the family Nudiviridae into four genera (Alphanudivirus to Deltanudivirus), as recently adopted by the International Committee on Taxonomy of Viruses (ICTV), but rather shows two main branches of the family that are further subdivided. Our data support a recent proposal to create two subfamilies within the family Nudiviridae, each subdivided into several genera.

Confirmation of Oryctes rhinoceros nudivirus infections in G-haplotype coconut rhinoceros beetles (Oryctes rhinoceros) from Palauan PCR-positive populations.

Coconut rhinoceros beetle (CRB), Oryctes rhinoceros, is a pest of palm trees in the Pacific. Recently, a remarkable degree of palm damage reported in Guam, Hawaii, Papua New Guinea and Solomon Islands has been associated with a particular haplotype (clade I), known as "CRB-G". In the Palau Archipelago, both CRB-G and another haplotype (clade IV) belonging to the CRB-S cluster coexist in the field. In this study, more than 75% of pheromone trap-captured adults of both haplotypes were Oryctes rhinoceros nudivirus (OrNV)-positive by PCR. No significant difference in OrNV prevalence between the haplotypes was detected. In PCR-positive CRB-G tissue specimens from Palau, viral particles were observed by electron microscopy. Hemocoel injection of CRB larvae with crude virus homogenates from these tissues resulted in viral infection and mortality. OrNV isolated from Palauan-sourced CRB was designated as OrNV-Palau1. Both OrNV-Palau1 and OrNV-X2B, a CRB biological control isolate released in the Pacific, were propagated using the FRI-AnCu-35 cell line for production of inoculum. However, the OrNV-Palau1 isolate exhibited lower viral production levels and longer larval survival times compared to OrNV-X2B in O. rhinoceros larvae. Full genome sequences of the OrNV-Palau1 and -X2B isolates were determined and found to be closely related to each other. Altogether these results suggest CRB adults in Palau are infected with a less virulent virus, which may affect the nature and extent of OrNV-induced pathology in Palauan populations of CRB.

Nudivirus Sequences Identified from the Southern and Western Corn Rootworms (Coleoptera: Chrysomelidae).

Analysis of pooled genomic short read sequence data revealed the presence of nudivirus-derived sequences from U.S. populations of both southern corn rootworm (SCR, Diabrotica undecimpunctata howardi Barber) and western corn rootworm (WCR, Diabrotica virgifera virgifera LeConte). A near complete nudivirus genome sequence was assembled from sequence data for an SCR population with relatively high viral titers. A total of 147,179 bp was assembled from five contigs that collectively encode 109 putative open reading frames (ORFs) including 20 nudivirus core genes. In contrast, genome sequence recovery was incomplete for a second nudivirus from WCR, although sequences derived from this virus were present in three geographically dispersed populations. Only 48,989 bp were assembled with 48 putative ORFs including 13 core genes, representing about 20% of a typical nudivirus genome. Phylogenetic analysis indicated that both corn rootworm nudiviruses grouped with the third known nudivirus of beetles, Oryctes rhinoceros nudivirus in the genus Alphanudivirus. On the basis of phylogenetic and additional analyses, we propose further taxonomic separation of nudiviruses within Alphanudivirus and Betanudivirus into two subfamilies and five genera. Identification of nudivirus-derived sequences from two species of corn rootworm highlights the diversity of viruses associated with these agricultural insect pests.

Chromosomal scale assembly of parasitic wasp genome reveals symbiotic virus colonization.

Endogenous viruses form an important proportion of eukaryote genomes and a source of novel functions. How large DNA viruses integrated into a genome evolve when they confer a benefit to their host, however, remains unknown. Bracoviruses are essential for the parasitism success of parasitoid wasps, into whose genomes they integrated ~103 million years ago. Here we show, from the assembly of a parasitoid wasp genome at a chromosomal scale, that bracovirus genes colonized all ten chromosomes of Cotesia congregata. Most form clusters of genes involved in particle production or parasitism success. Genomic comparison with another wasp, Microplitis demolitor, revealed that these clusters were already established ~53 mya and thus belong to remarkably stable genomic structures, the architectures of which are evolutionary constrained. Transcriptomic analyses highlight temporal synchronization of viral gene expression without resulting in immune gene induction, suggesting that no conflicts remain between ancient symbiotic partners when benefits to them converge.

Transcription Profile and Genomic Variations of Oryctes Rhinoceros Nudivirus in Coconut Rhinoceros Beetles.

Oryctes rhinoceros nudivirus (OrNV) is a double-stranded DNA (dsDNA) virus which has been used as a biocontrol agent to suppress the coconut rhinoceros beetle (Oryctes rhinoceros) in Southeast Asia and the Pacific Islands. A new wave of O. rhinoceros incursions in Oceania is thought to be related to the presence of low-virulence isolates of OrNV or virus-tolerant haplotypes of beetles. In this study, chronically infected beetles were collected from Philippines, Fiji, Papua New Guinea (PNG), and the Solomon Islands (SI). RNA sequencing (RNA-seq) was performed to investigate the global viral gene expression profiles and for comparative genomic analysis of structural variations. Maximum likelihood phylogenic analysis indicated that OrNV strains from the SI and Philippines are closely related, while OrNV strains from PNG and Fiji formed a distinct adjacent clade. We detected several polymorphic sites with a frequency higher than 35% in 892 positions of the viral genome. Nonsynonymous mutations were detected in several hypothetical proteins and 15 nudivirus core genes, such as gp034, lef-8, lef-4, and vp91 We found limited evidence of variation in viral gene expression among geographic populations. Only a few genes, such as gp01, gp022, and gp107, were differentially expressed among different strains. Additionally, small RNA sequencing from the SI population suggested that OrNV is targeted by the host RNA interference (RNAi) response with abundant 21-nucleotide small RNAs. Some of these genomic changes are specific to the geographic population and could be related to particular phenotypic characteristics of the strain, such as viral pathogenicity or transmissibility, and this requires further investigation.IMPORTANCE Oryctes rhinoceros nudivirus has been an effective biocontrol agent against the coconut rhinoceros beetle in Southeast Asia and the Pacific Islands for decades. The recent outbreak of these beetles in many South Pacific islands has had a significant impact on livelihoods in the region. It has been suggested that the resurgence and spread of the pest are related to the presence of low-virulence isolates of OrNV or virus-tolerant haplotypes of beetles. We examined viral genomic and transcriptional variations in chronically infected beetles from different geographical populations. A high number of polymorphic sites among several geographical strains of OrNV were identified, but potentially only a few of these variations in the genome are involved in functional changes and can potentially alter the typical function. These findings provide valuable resources for future studies to improve our understanding of the OrNV genetic variations in different geographic regions and their potential link to virus pathogenicity.

Visualizing Nudivirus Assembly and Egress.

Enveloped viruses hijack cellular membranes in order to provide the necessary material for virion assembly. In particular, viruses that replicate and assemble inside the nucleus have developed special approaches to modify the nuclear landscape for their advantage. We used electron microscopy to investigate cellular changes occurring during nudivirus infection and we characterized a unique mechanism for assembly, packaging, and transport of new virions across the nuclear membrane and through the cytoplasm. Our three-dimensional reconstructions describe the complex remodeling of the nuclear membrane necessary to release vesicle-associated viruses into the cytoplasm. This is the first report of nuclear morphological reconfigurations that occur during nudiviral infection.IMPORTANCE The dynamics of nuclear envelope has a critical role in multiple cellular processes. However, little is known regarding the structural changes occurring inside the nucleus or at the inner and outer nuclear membranes. For viruses assembling inside the nucleus, remodeling of the intranuclear membrane plays an important role in the process of virion assembly. Here, we monitored the changes associated with viral infection in the case of nudiviruses. Our data revealed dramatic membrane remodeling inside the nuclear compartment during infection with Oryctes rhinoceros nudivirus, an important biocontrol agent against coconut rhinoceros beetle, a devastating pest for coconut and oil palm trees. Based on these findings, we propose a model for nudivirus assembly in which nuclear packaging occurs in fully enveloped virions.

Endogenous viral elements integrated into the genome of the soybean aphid, Aphis glycines.

Sequence analysis of the genomic DNA isolated from four biotypes of the soybean aphid, Aphis glycines (AG), revealed that in addition to the commonly observed retrovirus-related retrotransposons, viral sequences derived from multiple RNA and DNA viruses have integrated into the genome. Notably, sequences of more than 60 nudiviral genes were identified from de novo assembled DNA contigs, and mapped to assembled genomic scaffolds of AG, indicating that an ancient nudivirus, named Aphis glycines endogenous nudivirus (AgENV), had integrated into the AG genome. Furthermore, sequences derived from a similar endogenous nudivirus, Melanaphis sacchari endogenous nudivirus (MsENV), were identified from the genomic scaffolds of the sugarcane aphid, Melanaphis sacchari. Analysis of transcriptome and small RNA sequence data derived from AG did not provide evidence for transcription of the integrated AgENV genes. Hence, the genes of AgENV may be present as pseudogenes. Phylogenetic analysis based on nudivirus core genes indicated that these aphid ENVs belong to the genus Alphanudivirus.

Nudivirus Remnants in the Genomes of Arthropods.

Endogenous viral elements (EVEs), derived from all major types of viruses, have been discovered in many eukaryotic genomes, representing "fossil records" of past viral infections. The endogenization of nudiviruses has been reported in several insects, leading to the question of whether genomic integration is a common phenomenon for these viruses. In this study, genomic assemblies of insects and other arthropods were analyzed to identify endogenous sequences related to Nudiviridae. A total of 359 nudivirus-like genes were identified in 43 species belonging to different groups; however, none of these genes were detected in the known hosts of nudiviruses. A large proportion of the putative EVEs identified in this study encode intact open reading frames or are transcribed as mRNAs, suggesting that they result from recent endogenization of nudiviruses. Phylogenetic analyses of the identified EVEs and inspections of their flanking regions indicated that integration of nudiviruses has occurred recurrently during the evolution of arthropods. This is the first report of a comprehensive screening for nudivirus-derived EVEs in arthropod genomes. The results of this study demonstrated that a large variety of arthropods, especially hemipteran and hymenopteran insects, have previously been or are still infected by nudiviruses. These findings have greatly extended the host range of Nudiviridae and provide new insights into viral diversity, evolution, and host-virus interactions.

Chalcid wasp paleoviruses bridge the evolutionary gap between bracoviruses and nudiviruses.

Polydnaviruses are obligate mutualists of parasitoid wasps and are divided into two genera, Bracovirus and Ichnovirus. Bracoviruses are thought to originate from a single integration of an ancestral nudivirus into the ancestor of microgastroid complex ~100 million years ago. However, all the known nudiviruses are only distantly related to bracoviruses, and much remains obscure about the origin of bracoviruses. Here we employ a paleovirological method to screen endogenous nudivirus-like elements across arthropods. Interestingly, we identify many endogenous nudivirus-like elements within the genome of Eurytoma brunniventris, a species of the Chalcidoidea superfamily. Among them, we find 14 core gene sequences are likely to be derived from a betanudivirus (designated EbrENV-ß), suggesting that betanudivirus has been circulating in parasitoid wasps. Phylogenomic analysis suggests that EbrENV-ß is the known closest relative of bracoviruses. Synteny analyses show the order of core genes is not well conserved between EbrENV-ß and nudiviruses, revealing the dynamic nature of the evolution of nudivirus genome structures. Our findings narrow down the evolutionary gap between bracoviruses and nudiviruses and provide novel insights into the origin and evolution of polydnaviruses.